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ABSTRACT 


Naval operational weather products, especially EM prediction, often fail to 
capture key meteorological features in the boundary layer due to the absence of 
high-resolution profile data. Quad-rotor small unmanned aerial systems (sUAS) 
with meteorological sensing capabilities may provide a solution for acquiring 
vertical profile data at sea. With the use of an RS92 radiosonde, InstantEye 
sUAS, and a calibrated sensing tower, we analyzed the effects a quad-rotor had 
on the accuracy of temperature and pressure profiles in the surface layer. 

In unstable atmospheres, temperature measurements made in the surface 
layer are as accurate as the manufacturer’s claimed accuracy. In stable 
atmospheres, mixing occurs below 1.3 m, while above 2 m, sampled air comes 
from as much as 2 m aloft. The InstantEye’s rotors contribute to the variation in 
temperature measurement, and this effect is strongest near the surface. The 
variations introduced by the prop-wash helped suppress natural variations from 
turbulent fluctuations, but are still present above the surface layer. The 
InstantEye also introduces a 0.12 hPa pressure bias while in flight due to rotor- 
induced lift. 
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I. INTRODUCTION 


A. UNMANNED AERIAL SYSTEMS 

The use of unmanned aircraft in military appiications is nothing new. 
Remoteiy piloted airplanes can be traced back to as early as World War I in 
which British researcher A. M. Low developed his “Aerial Target” to fly ordinance¬ 
laden planes into enemy positions and aircraft (Tayior 1977). While the concept 
appealed to military leaders, the technology was severely lacking to provide any 
decisive advantage on the battiefieid. Years later, in Worid War II, more robust 
aircraft with superior controi systems saw action in combat, but still failed to 
demonstrate their necessity as part of a military arsenal (Keane and Carr 2013). 

Despite these early failures, interest in unmanned aircraft persisted in 
miiitary circles throughout the years. Today advances in communications, signai 
processing, and manufacturing techniques have enabled aircraft designers to 
provide miiitaries with effective mission-capabie remotely-piioted assets. 
Common names for these include drones, unmanned aerial vehicles (UAV), 
remotely piloted aircraft (RPA), unmanned aeriai system (UAS), or small LIAS 
(sUAS). A UAS is defined as, “a system whose components inciude the 
necessary equipment, network, and personnei to control an unmanned aircraft” 
(JCS 2010). For the purposes of this document, we will use only the term UAS 
and sUAS to describe unmanned aircraft, as the definition provides a broader 
perspective for describing all components related to unmanned flight. 

Today the U.S. miiitary maintains the iargest investment in UAS 
technoiogy compared to any other organization. As of Juiy 1, 2013, the DOD 
maintained nearly 11,000 UAS in its inventory, allocating $3.8 biiiion from the 
FY2014 budget. Despite projected budgetary cutbacks, funding is expected to 
grow by 17% through FY2018 (Pentagon 2014). Such commitment to UAS 
integration demonstrates the utiiity UAS provide on the battiefieid. The platforms 
in use today perform a variety of missions including intelligence, surveillance, 
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and reconnaissance (ISR), as well as combat strike missions. The focus of future 
UAS military technology will be primarily centered on ensuring dominance in 
contested airspace, and superiority in anti-access area denial (A2/AD) scenarios. 
This not only requires improvements to communication and sensing networks for 
tactical efficiency, but also necessitates an improvement in environmental 
awareness to better assess risk and exploit certain weather phenomena such as 
electromagnetic (EM) propagation features. 

B. WEATHER AND EM PROPAGATION PREDICTION 

Understanding the weather is a critical component in combat 
effectiveness. Every military asset can be adversely affected by unfavorable 
weather conditions; therefore, operational commanders depend on accurate 
forecasts. This is particularly challenging for the Navy where operations take 
place over open water and routine meteorological measurements are sparse. In 
the past, naval forecasters, referred to as Aerographer’s Mates (AGs), utilized 
weather balloons to sample the vertical structure of the atmosphere. This not 
only assisted in developing forecasts for daily operations, but also improved 
shore-based computer weather prediction models since the data could easily be 
forwarded and integrated into the model. Unfortunately, the Navy’s operational 
weather balloon program was terminated in 2011 (OPNAV N2/N6F5 2010), 
causing tremendous uncertainty in vertical structure of the atmosphere in crucial 
regions. Since then AGs have had to rely exclusively on model data, satellite 
images, and surface observations to perform their duties while at sea. 

A consequence of the knowledge gap in vertical atmospheric structure is 
the inability to create accurate EM predictions products. These products predict 
how various communication and radar systems will perform based on given 
atmospheric conditions. The atmosphere, particularly near the surface, can 
experience large diurnal and daily changes. These variations influence the 
effective range over which a radar system can detect a target. An aircraft might 
have an extended communication range one day, and a severely limited range 
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the next. EM signals also have the potential of being intercepted or jammed, so 
understanding the performance parameters of systems belonging to both friendly 
and adversarial forces is crucial in wartime. Recognizing the importance of EM 
propagation in future warfare, the Chief of Naval Operations ADM Jonathan 
Greened stated in 2011, “Control of information, much of it through the 
electromagnetic spectrum, is already growing more important than the control of 
territory in modern warfare” (OPNAV N2/N6 2012). EM waves are the means 
with which information is exchanged on the battlefield, and knowledge of how 
they propagate can mean the difference of exploiting an adversary’s activities, or 
being exploited by them. 

EM signals In radio frequencies are influenced by three atmospheric 
properties: temperature, pressure, and relative humidity (Ko et al. 1983). 
Variations of these parameters are used by the military’s Advanced Refractive 
Effects Prediction System (AREPS) to analyze how various radar and 
communication systems are expected to operate (SPAWAR 2006). While 
AREPS has the capability for the user to manually input atmospheric parameters 
(e.g., from a weather balloon), the atmospheric variables are now usually 
obtained through the Coupled Ocean/Atmosphere Mesoscale Prediction System 
(COAMPS™) weather model. However, poor vertical resolution in COAMPS™ 
output presents a major disadvantage. Many assumptions about the vertical 
profile are made, which can have a dramatic effect on EM propagation products. 
In the atmospheric boundary layer (ABL), or the layer approximately 100 to 3000 
meters above ground (Stull, 1998), surface fluxes can dramatically affect 
temperature and humidity profiles, often missed by COAMPS™. Consequently, 
AREPS may fail to accurately capture significant propagation features. 

C. WEATHER SENSING WITH UAS 

With the increasing use of military UAS, it seems logical to employ these 
aircraft to sample the atmosphere. Since most UAS are smaller and lighter than 
traditional aircraft due to lack of a human pilot and the accompanying safety 
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systems, they tend to be more vulnerable to adverse weather conditions such as 
wind gusts and icing. Atmospheric sensors installed onboard could not only serve 
to alert operators of dangerous conditions, but the data could also be collected 
for weather model inputs to improve future forecasts. 

Some UAS, such as the NOAA-owned Global Hawk, are used exclusively 
for high altitude atmospheric measurements (Braun 2013). These aircraft have 
been outfitted with an array of meteorological sensors used to conduct research 
in tropical cyclones. The DOD is also mounting meteorological sensors, and 
several noteworthy experiments have occurred in recent years. BAE Systems’ 
Manta and Boeing’s Scan Eagle were tested in 2011 and 2012, demonstrating 
the first known UAS measurements of heat, water vapor and momentum fluxes 
at low altitudes (Reineman 2013). In another experiment as part of exercise 
Trident Warrior 2013 (TW13), the Scan Eagle was again tested to capture 
boundary layer profiles, and successfully demonstrated that atmospheric 
measurements collected by UAS could be ingested into COAMPS improving 
short-range forecasts (NRL 2015). While this may at first seem promising, TW13 
also revealed that due to Scan Eagle’s safety requirements restricting low-level 
flight, significant EM propagation features were not captured when compared to a 
tethered sensor in the lowest 100 meters. 

UAS such as the Manta or Scan Eagle are fairly large enough fixed-wing 
aircraft to require sophisticated launch and recovery apparatuses and a team of 
skilled workers to operate and maintain them. Their procurement cost is also 
substantial, thus explaining any apprehension toward near-surface flight. A 
commander in possession of one or more UAS must also balance a set of 
demanding mission requirements with all available resources. For these reasons, 
a fixed-wing UAS already in service is not an ideal candidate to be used as a 
dedicated meteorological sensor. sUAS, on the other hand, may be well suited to 
serve in this capacity, and may be a feasible alternative to their larger and more 
expensive counterparts. 
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Small quad-rotor helicopters, also known as quad-copters, are a class of 
sUAS. Quad-rotor sUAS, as their name implies, have four rotors to generate lift, 
making them very stable aircraft. They can also be very compact, measuring less 
than one-half meter in width and weighing only a couple kilograms. Unlike their 
larger fixed-wing counterparts, a single person with very limited training can 
operate quad-rotor sUAS. Due to large-scale commercialization of quad-rotors 
over the past decade, their costs have significantly declined to only hundreds or 
a couple thousand dollars; a mere fraction of the cost associated with traditional 
DOD UAS. The simplicity and affordability of quad-rotor sUAS make it feasible to 
employ them on a large scale as dedicated meteorological sensing devices. A 
vast network of worldwide sensors in the ABL could significantly improve model 
products and operational EM propagation forecasting techniques. 

D. RESEARCH OBJECTIVES 

The use of small quad-rotor sUAS for meteorological sampling is a fairly 
novel concept with very little literature existing on the subject. The aim of this 
document is to investigate how the performance of meteorological sensors 
attached to a miniature quad-rotor sUAS is affected while in flight. Specifically, 
temperature and pressure datasets were collected during Joint Interagency Field 
Exercises (JIFX) in 2014 for both stable and unstable atmospheres. Our analysis 
is an investigation of sensing performance under these two conditions. 
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II. BACKGROUND AND CONSIDERATIONS 


The aerodynamics and various motions of air surrounding an aircraft with 
a spinning rotor are complex. Four rotors operating in an equilateral array add 
another layer of complexity to the problem. Furthermore, atmospheric conditions 
in the boundary layer are also turbulent and difficult to quantify from single point 
measurements. Interaction between air and the earth’s surface drive many 
processes in fluid motions and gradient distributions. To understand how a quad- 
rotor sUAS in flight affects the ABL meteorological conditions that it is supposed 
to measure, it is important to first develop some background information 
pertaining to these two subjects. 

A. QUAD-ROTOR MECHANICS 

The principles governing quad-rotor performance are related to the 
helicopter. The following discussion on rotor flight performance is largely based 
on the work of Leishman (2006). 

1. Momentum Theory 

A general helicopter undergoes many different flight configurations such 
as hovering, forward motion, and climb and descent. The mass flow through the 
rotor is a complicated process of varying velocities accompanied by a sequence 
of pressure-induced vortices in the slipstream (or wake) that differ between each 
configuration of flight. However, rotor performance can be analyzed through a 
simplified method in which assumptions are made about the model. This method 
is well adopted in the field of aerodynamics known as the Rankine-Froude 
momentum theory. 

The model assumes that the air is incompressible and inviscid, which 
generates no shear within the fluid. The flow through the rotor is one¬ 
dimensional, steady and uniform, and does not impart a rotational component to 
the traversing fluid. While these assumptions will inevitably introduce errors into 
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the results, the simplification outweighs these errors, which can be mitigated 
through an appropriate correction factor applied to the final result (Layton 1984). 
Momentum theory also assumes that the rotor is an infinitesimally thin actuator 
disk, which supports the thrust generated by the rotational motion of the blades. 
Work is done on the rotor from torque applied by the rotor shaft, which is 
converted into a net gain of kinetic energy within the rotor wake. 

a. Stationary Flight 

A hovering helicopter is the most basic configuration to analyze 
performance. There is zero vertical or horizontal velocity, and the rotor flow field 
is said to be azimuthally axisymmetric. Consider a finite control volume S 
containing the rotor and surrounding air volume as shown in Figure 1. The top 
plane is at a distance far above the rotor and the fluid is quiescent (i.e., = 0). 

Plane 1 and plane 2 represent the area just above and below the rotor, 
respectively. Given an infinitesimally thin rotor plane, the fluid velocity is equal 
through these two planes at V^ + v.', however, a pressure differential exists 

across the rotor plane enabling the thrust needed to counter the weight of the 
aircraft. Far below the rotor plane is the wake plane, denoted as oo, with the 

wake velocity V^^ + w. From the unit normal area vector, denoted by dS, and the 

conservation of mass we can express the mass flux through the control volume 
as 


jj^pdVdS = 0, ( 1 ) 

where p is air density. 

Put another way, the mass flowing into the control volume is equal to the 
mass flowing out. The mass flow rate m must also be consistent through any 
portion of the control volume. Hence, 

m = jj pV dS = jj pV dS = jj^pV dS, (2) 
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which is further reduced by the one-dimensional incompressible flow assumption 


m = pA^w = pAVf. 


( 3 ) 





/ 


\ t / 

I Weight j 


00 


Vc+ w 

Figure 1. Flow model for momentum theory analysis of a rotor in 
hovering flight (from Leishman 2006). 

The work done on the rotor is equal to the kinetic energy imparted into the 
slipstream per unit time, also known as power Tv^, where T is the upward force 

or thrust exerted by the rotor. Power can then be expressed as 



( 4 ) 


However, the velocity at plane 0 is nominally zero, eliminating the right hand 
term. Thus the equation reduces to 


Tv. — —mw 


2 


( 5 ) 
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V. = —w. 
2 


(6) 


Combining Equations (3) and (6) provides the relationship 

pAVi = pA^w = pA^ (2v,.). (7) 

Stated another way, the wake velocity area is half the area of the rotor disk. 

From the velocity-area relationship, fluid motion about a rotor plane can be 
approximately depicted. However, with three additional rotors on a quad-rotor 
sUAS the flow becomes a bit more complicated. Applying the assumptions of 
momentum theory to a quad-rotor during hovering flight, a graphical 
representation of fluid flow can be constructed as shown in Figure 2. Given a 
wake velocity of approximately 4 m/s, and four rotors placed in a square pattern 
each with a radius of 7.5 cm, the control volume is now unique. There is a single 
quiescent plane at y = 0 diverging into four separate slipstreams as it passes 
through the rotor plane as depicted by the green square plane of Figure 2. To 
minimize slipstream influence on an attached meteorological sensor, ideal 
placement would be somewhere below the rotor plane and outside of the 
slipstreams. 
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Graph ica] Depict ion of Wind Velocity Through a Ouadrotor 



Figure 2. Graphical depiction of fluid velocity through a quad-rotor sUAS 

b. Climb and Descent 

While hovering flight is relatively straightforward to analyze, more complex 
configurations need to be considered for a more comprehensive understanding 
of fluid flow through a rotor. Climbing and descending are necessary maneuvers 
for any aircraft, especially one designed to profile the lower atmosphere. When a 
rotor experiences vertical motion, the air far above and below can also be seen 
as moving vertically as experienced from the rotor’s point of view. This addition of 
velocity changes the various forces involved in flight. 

In a climbing configuration ^^>0, and there is a net upward force. The 
mass flow rate can be rewritten as 


m = pA^ (Vp -I- w) = pA{V(, + Vi) (8) 

Likewise, the thrust produced Is now 


T = m(y^ +w)- mV^ = mw 
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( 9 ) 





































( 10 ) 


T = pA{Vc + v,.)w = 2pA{Vc + v,.)v,.. 

Now let vl equal one-half the thrust over density and area, 


2pA 


= vl=vt+V,v,_ 


( 11 ) 


Rewriting (11) as a quadratic, it has the solution 


Vjv.^ 

o = l^| +-1^1-1 


( 12 ) 




(13) 


From our flow model, v, /v;,must always be positive for climbing flight, 
eliminating the negative solution. The increase in climb velocity implies a 
decrease in the rotor-induced velocity, which is referred to as the normal working 
state of the rotor (Figure 3). The limit is reached, however, at hovering flight just 
before descending. When -2v^<V^<0 the velocity at any plane can be either 

upward or downward creating a complex and turbulent flow pattern. Momentum 
theory is no longer applicable in this region. The slipstream will exist entirely 
above the rotor flowing in reverse direction when V^<-2v^. The rotor now 
extracts power from the slipstream and is referred to as the windmill brake state. 
Since during the descent. Equations (8) through (13) can be reassessed 

to form the solution 


V,. 




(14) 


Again, there are two solutions, but only the negative version is valid to 
ensure v, /v;, <1. This is logical since a reduction in induced rotor velocity results 

in a faster descent. Theoretically, rotor-induced velocity contributions should be 
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negligible to ensure the accuracy of atmospheric measurements. Favorable 
conditions exist during periods of climbing and of rapid decent. Slow descent is 
unfavorable since the assumptions of momentum theory are no longer valid, and 
the stream velocity is unpredictable, creating turbulent conditions within the fluid. 



Climb velocity ratio, 


Figure 3. Induced velocity as a function of climb and descent from 
momentum theory (from Leishman 2006). 


c. Horizontal Motion 

In the previous sections on helicopter dynamics, we were able to merge 
momentum theory from a traditional rotorcraft to a quad-rotor model. However, 
the horizontal motion of a quad-rotor is fundamentally different than a traditional 
helicopter system. Traditional helicopters only have two rotors and are used in 
one of two ways, (1) both rotors are used to generate thrust, but rotate in 
opposite directions to counter the rotor-induced torque, or (2) one rotor is used to 
generate thrust, and a smaller rotor is used to counter the torque acting on the 
fuselage. Horizontal motion is accomplished by changing the angle of attack 
(AoA) of the rotor plane and controlled through the adjustment of the rotor blade 
pitch and angle relative to the drive shaft. 
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By contrast, miniature quad-rotor sUAS have rigid rotor blades that do no 
vary in pitch and angle. The four rotors are equally separated in a square pattern, 
and each spins in the opposite direction as the rotor adjacent to it, balancing the 
four torque forces (Figure 4). There are only four control inputs (i.e., the rotation 
speed of each rotor) to regulate the aircraft’s pitch, roll and yaw. 



\ 


e3 


e2 


Figure 4. Forces acting on a quad-rotor aircraft 
(from Tayebi and McGilvray 2006). 


Horizontal movement is achieved by changing the pitch or roll of the 
airframe by changing the speed of one or more rotors. The thrust generated by 
the four rotors is no longer perpendicular to the ground since the rotor plane is 
tilted at an angle. As a consequence, the flow through the rotors is no longer 
axisymmetric since the thrust generated is partially transferred to into a 
horizontal component for movement. The streamline in this configuration is also 
non-perpendicular to the ground and is angled in the opposite direction of 
movement. The angle of the streamline could affect atmospheric measurements 
if the sensor should find its way into the rotor-induced velocity field. Therefore, a 
meteorological sensor would likely perform best when positioned near the front of 
the aircraft. 
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2 . 


Vortices 


Vortices are a byproduct of lift generated by an airfoil. The pressure 
differential that exists on either side of the airfoil creates lift. However, at the wing 
tip the solid boundary disappears, and the difference of pressure must be 
equalized. The high pressure under the airfoil flows toward the lower pressure 
above the wing, and the result is a wake vortex. A rotor is an airfoil rotating at 
high speeds, and as such Leishman explains that a set of vortices are created 
that follow the rotation of movement. Since the slipstream velocity is downward in 
a hover or climbing configuration, the vortices will also drift downward creating 
turbulence in their path (Figure 5). 


Rotor thrust 



Tip vortices convect 
along slipstream boundary 


Figure 5. Velocity field of a two-bladed rotor operating in a hover 

(from Leishman 2006) 

Turbulence from vortices complicates the flow structure around the 
aircraft. This is troubling from a meteorological perspective as the rotational 
energy can promote mixing of the surrounding fluid affecting the accuracy of 
measurements. This topic is not well understood, and the extent to which mixing 
affects the environment is in part a goal of this experiment. 
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3 . 


Ground Effect 


Another phenomenon affecting helicopter performance is ground effect. 
This occurs when the aircraft operates at very low altitudes, and the earth’s 
surface constrains the development of its rotor wake. This has several 
implications for low-level flight, but predominantly affects the power-to-thrust 
ratio. When ground effect is present, less power is needed to achieve the 
necessary thrust required to sustain flight (Leishman 2006). Shown in Figure 6, 
the ground is a boundary to the streamline flow in the presence of ground effect. 
As a result, the slipstream rapidly expands as it approaches the ground, altering 
the slipstream and rotor-induced velocities. From momentum theory, these 
velocities are related to the power and thrust requirements. 


(a) Out of ground effect (OGE) 


(b) In ground effect (IGE) 


No effect on 
wake near rotor 


Some effect 
on wake near ^ 
ground 


4 No effect on rotor 
thrust or power 



Large effect on rotor thrust 
and/or power 

Rotor height ^ ^ 

off the ground 

Rotor 


Ground 



Figure 6. Wake behavior: (a) Out of ground effect, (b) In ground effect 

(from Leishman 2006) 


Ground effect’s influence on rotor performance varies by altitude. A rotor 
closer to the ground will experience more slipstream restriction and have a higher 
thrust-to-power ratio than a rotor at a higher altitude. According to Leishman 
(2006), empirical observations suggest that significant performance effects due 
to ground effect occur below an altitude of one rotor diameter. To sufficiently 
minimize ground effect, an altitude of at least two rotor diameters is 
recommended. 
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B. SURFACE LAYER METEOROLOGY 


The atmosphere is composed of several horizontal layers. The region 
where weather occurs is the troposphere, which extends from the surface to an 
average altitude of 11 km. In the lowest portion of the troposphere is the 
boundary layer, defined by Stuii (1988) as, “That part of the troposphere that is 
infiuenced by the presence of the earth’s surface, and responds to surface 
forcings with a timescaie of about an hour or iess.” The boundary layer’s 
tendency to change rapidly has significant impacts to iow-levei meteoroiogy. The 
ground absorbs much of the incoming soiar radiation, and in turn the warming 
ground heats the atmosphere above it through turbulent fluxes and emitting 
radiation. The energy transfer between the earth’s surface and iower atmosphere 
is distributed through a process called turbulent transport in which warm eddies 
develop near the ground and are lifted aioft by their buoyancy. 

The lowest part of the boundary layer is called the surface layer. Stull 
(1988) refers to this layer as, “The region at the bottom of the boundary layer 
where turbulent fluxes and stress vary by less than 10%, decreasing above that 
in the ABL.” In other words, the energy transfer is strongest very close to the 
surface. The surface layer is important to understand since the majority of the 
data collected in the experiment was sampled in this region. 


1. Potential Temperature 

Potential temperature (0) is a conserved variable in an adiabatic process, 
accounting for changes in pressure (Petty 2008). When comparing a sample of 
vertical temperatures it is important to use potentiai temperature since according 
to the ideal gas law temperature and pressure are interdependent. The 
temperature can be corrected for the pressure change by considering potentiai 
temperature defined as 


e = T 


rp\0.2S6 


KPJ 


(15) 
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where T is temperature in degrees Kelvin, P is the air pressure in hPa, and is 

a reference pressure normally 1000 hPa. However, boundary layer work requires 
using the surface pressure so (15) can be approximated by 


e = T + 



z, 


(16) 


where z is the height in meters above the ground, and (g/C^) is the dry 
adiabatic lapse rate of 0.0098 K/m, composed of g, the acceleration due to 
gravity, and C^, heat capacity at constant pressure (Stull 1988). 


2. Humidity 

Humidity is an important quantity affecting EM propagation through the 
atmosphere. The amount of moisture a parcel of air can hold is directly 
proportional to that parcel’s temperature. Saturation will occur If the moisture 
content of a parcel equals its temperature limit. Often relative humidity (RH) is 
used to describe this relationship as a percentage, and is calculated 


RH = 



(17) 


where e is the actual vapor pressure, and e^{T) is the saturation vapor pressure. 

Relative humidity is commonly measured and reported by weather sensors. 
However, when comparing a set of vertical measurements In the presence of a 
temperature gradient, specific humidity (^, the ratio of mass of water to the mass 
of air) is more convenient to use since it is a conserved quantity. Specific 
humidity is calculated 


w 

q = -— > 

1 +w 


(18) 


where w is the mixing ratio, defined as the ratio between the mass of water 
vapor to the mass of dry air. 
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III. EXPERIMENT 


A. METHODOLOGY 

Testing was performed exclusively at McMillan Airfield, located at Camp 
Roberts, CA, as part of the Joint Interagency Field Experimentation (JIFX) 
program in May, August, and November 2014. An InstantEye Mark 2 Gen 3, 
produced by Physical Sciences Inc. (PSI), was flown in an open field next to a 
calibrated meteorological observation tower. The tower was in constant operation 
over the course of each multi-day JIFX measuring and recording pressure, wind, 
soil parameters, and seven separate levels of temperature and relative humidity. 
The tower was strategically placed away from obstructions in an open field to 
minimize localized weather effects. 

To measure atmospheric parameters onboard the sUAS, a Vaisala RS92 
radiosonde was attached to the underside of the InstantEye. Originally designed 
as an instrument for weather balloons, the radiosonde proved to be an excellent 
instrument for this experiment due to its light weight, comprehensive sensing 
package, and communication system. When in operation, the radiosonde 
samples temperature, pressure, and relative humidity while transmitting its data 
in real time to a computer located in a nearby hangar. 

Test flights were conducted at various times to collect a diverse data set. 
Of the flights that were later analyzed, five occurred during daylight hours when 
the sun’s warming effect produced an unstable temperature profile near the 
observation tower, and two were conducted just before sunrise when the effect of 
radiative cooling was at its maximum, producing a very stable temperature 
profile. Security measures in place at Camp Roberts made access to the tower 
before 6:00 AM difficult, limiting the amount of testing that could be conducted in 
a stable atmosphere. 
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Figure 7. Physical Sciences Inc. InstantEye with 
Vaisala Radiosonde Attached 


) 

■# 



Figure 8. Meteorological Observation Tower with 
seven sensor levels and one wind level. 
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A typical test flight involved hovering the aircraft at a level corresponding 
to one of the seven sensors on the observation tower. While the radiosonde 
collected atmospheric information in flight, the corresponding tower sensor 
served as a reference instrument recording temperature and relative humidity 
unaffected by processes produced by the aircraft. Synchronization between the 
tower and radiosonde was crucial for proper analysis, and achieved through the 
use of an audio recorder. Recordings were used to document time offsets 
between the observation tower and the radiosonde computer, as well as provide 
exact timing of aircraft vertical position. 

B. INSTANTEYE MK 2 GEN 3 

PSI’s InstantEye was designed for tactical level ISR operations in urban 
environments giving small ground units situational awareness through the use of 
its onboard cameras. This sUAS was chosen for the experiment due to its 
already existing employment within the DOD, active support from PSI. It is likely 
that, any quad-rotor sUAS capable of carrying the radiosonde payload could also 
be used in a similar manner. 

The InstantEye measures 25 cm across between rotor axes, with each of 
its four rotors measuring 20.3 cm diameter, and weighs only 320 g without 
battery or payload. It is capable of operating in all weather, temperatures ranging 
from -18 to 49 C, and can endure winds of 13 m/s. The endurance rating is 
dependent on payload weight, battery charge, and airfoil efficiency, but typical 
flights during the experiment lasted about 15 to 20 minutes before battery 
replacement was needed. 

Built into the InstantEye are three cameras at different angles, GPS, and 
pressure sensor, which send information to the operator via the flight controller 
screen. The camera image is displayed in real time on the controller serving as a 
tool for navigation when visual contact is lost. The GPS device provides real time 
location coordinates of the aircraft, and also serves to hold the aircraft’s lateral 
position under autonomous flight making it less susceptible to unwanted 
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movement due to prevailing and gusting winds. The GPS information is part of 
the InstantEye’s internal circuitry and was not recorded during the experiment. 
The pressure data measured from the pressure sensor is unavailable to the user, 
but the information is converted to altitude and displayed in feet. The InstantEye 
responds to slight changes in pressure to control its climb rate and maintain its 
altitude. Like the GPS data, the InstantEye pressure measurements were not 
recorded during the experiment and only used for avionics control. 




Figure 9. PSI InstantEye Front and Bottom View 
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c. 


RADIOSONDE 


For the purpose of the experiment, Vaisala’s RS92 radiosonde was an 
ideal meteorological sensing device to be used onboard the aircraft. Vaisaia has 
been a trusted name for reiiable and accurate meteoroiogical sensors for 
decades, and the RS92 is a model is used woridwide for verticai atmospheric 
sampling attached to weather balloons. According to Vaisaia’s datasheet, the 
sensing package inciudes a single wire capacitance temperature sensor, duai 
heated thin-fiimed capacitance relative humidity sensors, and a micromechanicai 
silicon pressure sensor. Capacitive sensors are not infiuenced by magnetic 
variations (Kruis et ai. 1998), and this is important since the sUAS is controlled 
by electric signais, which generate magnetic disturbances very close to the RS92 
sensors. 

The RS92 also contains a GPS device to track iaterai movements aloft 
used to calculate wind speed and direction whiie attached to a drifting bailoon. 
Since the sensor was attached to a controiied aircraft designed to resist wind 
variations, wind data coiiected from the RS92 during the experiment were 
irrelevant. It was thought that the GPS coordinates collected might be able to 
help determine iaterai movements of the aircraft during the experiment, but the 
resolution was too coarse to extract any meaningful information. 

Another feature of the RS92 is a fully digital transmission signal on the 
400-405 MHz meteorological band. In general, there are two methods of 
collecting data from a sensor attached to the aircraft, 1) the data is stored on a 
locai memory device, then extracted to a computer once the aircraft is recovered, 
or 2) the data is streamed wireiessly while in fiight to a recording device on the 
ground. Locally stored data as in option (1) presents some risk. If the aircraft 
crashes (as it did many times during the experiment), or is unrecoverable, there 
is the potential of also losing the data coiiected. For the experiment a iocaiiy 
stored device might be have been adequate since most of the flights took place 
over land at low levels, but in a real at-sea application aititudes would be much 

higher over open ocean furthering the risk of non-recoverabiiity. While a 
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specifically designed sensor meeting all the requirements for mission deployment 
is preferred, the RS92 is an affordable option and could be used immediately for 
fleet deployment, making it an ideal test instrument for this experiment. 



Figure 10. Vaisala Radiosonde RS92-SGP 


SENSOR 

OPERATING 

RANGE 

ACCURACY 

Temperature 

-90 to 60 C 

±0.5C 

Pressure 

3 to 1080 hPa 

± 1 hPa 

Relative Humidity 

Oto 100% 

±5% 


Table 1. Vaisala Radiosonde RS92-SGP specifications 

(from www.vaisala.com) 
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D. OBSERVATION TOWER 

The tower was constructed using a metallic tripod extending upward to 
6.2m where a Vaisala WXT520 weather station was mounted. For the purpose of 
the experiment only barometric pressure, wind speed, and wind direction were of 
any value by the WXT520. Separate temperature and relative humidity probes 
were installed along the vertical axis of the tripod with approximately logarithmic 
spacing. The spacing ensured a near linear temperature sampling distribution 
due to the rapid temperature differential caused by surface heat flux (Table 1). 


Tower Level 

Height (m) 

WXT520 

6.20 

7 

4.93 

6 

2.44 

5 

1.30 

4 

0.66 

3 

0.34 

2 

0.19 

1 

0.11 


Table 2. Observation tower sensor heights 


Seven Rotronic MP100H relative humidity sensors with accompanying 
Pt100 1/10 DIN RTD temperature sensors were encased in PVC piping attached 
to the tripod. The portion of housing wall enclosing the sensors was double 
layered, fitted with a shield to minimize contamination from solar heating. An 
aspirator was used in the bottom six levels ensuring a uniform 3.5 m/s flow over 
the sensors. The seventh level was a later addition after the tower’s original 
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construction, and because of this, a separate aspirator and ventilation path with 
an unknown flow rate exist in the highest temperature and relative humidity 
sensors. 

During the experiment an additional temperature probe with a faster 
response time was required to capture turbulent temperature variations. A 
Campbell Scientific FW05 fast response thermocouple was installed at level four 
outside the PVC housing exposed to the elements. According to FW05 brochure, 
solar contamination is negligible due to the small diameter of the FW05 wire 
(0.0127 mm), and a radiation shield is not required. Each of the tower sensors 
were linked to a Campbell Scientific CR3000 Micrologger, which recorded at a 
sample rate of 10 seconds (0.1 FIz). The accuracy of each instrument is 
summarized in Table 3. 


SENSOR 

MFR. 

MODEL 

OPERATING 

RANGE 

ACCURACY 

Pressure 

Vaisala 

WXT520 

600 tollOOhPa 

±0.5 hPa 
(0-30 C) 

Wind Speed 

Vaisala 

WXT520 

0 to 60 m/s 

±3% 

(10 m/s) 

Wind Direction 

Vaisala 

WXT520 

Oto 360° 

o 

CO 

+1 

Temperature (RTD) 

Rotronic 

PtIOO 

1/10 DIN 

-40 to 60 C 

± (0.3 -H 0.005*t) 

Temperature (TC) 

Campbell 

Scientific 

FW05 

-270 to 1000 C 

±1.7C 

Relative Humidity 

Rotronic 

MP100H 

Oto 100% 

± 1% 


Table 3. Observation tower sensor specifications 
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IV. ERROR ANALYSIS 


To determine whether the InstantEye prop-wash or other contamination 
effects have any influence over the accuracy of the onboard radiosonde sensors, 
a method must first be developed to determine an acceptable range of accuracy. 
No matter how carefully an experiment is carried out, every measurement 
contains within it some degree of error. This chapter is aimed at examining the 
various sources of errors and uncertainties introduced in the experiment, and to 
determine an appropriate correction and tolerance for measured values. 

A. RANDOM AND SYSTEMATIC ERRORS 

Experimental uncertainties exist in two basic categories, random and 
systematic. The former can be revealed through repeated measurements (e.g., 
averaging a number of independent measurements to approximate the true value 
being measured). However, systematic errors cannot be resolved through 
repeated measurements. These errors always exhibit some degree of uncertainty 
biased in one direction away from the true value being measured (e.g., a clock 
running consistently slow will always underestimate the true time in an 
experiment no matter how many trials are conducted) (Taylor, 1997). 

Almost all measurements are subject to both systematic and random 
errors, but it is often possible to correct for both types of errors to reduce the final 
uncertainties. Random errors are identified through repeated measurements, and 
statistical methods exist to estimate how much the uncertainty is reduced with 
more sampling. A series of repeated measurements can be averaged, and this 
value (J) is assumed to be more reliable than any individual measurement alone 
for determining the mean value. The standard deviation {a^) is a measure of the 

variance individual measurements deviate from the average and is a measure of 
the uncertainty of a single measurement. One standard deviation represents a 
68% confidence interval that the true value of x lies somewhere between T + 

However, if these measurements are truly random and independent from each 
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other, the total uncertainty can be reduced by calculating the standard deviation 
of the mean (cTj), and is given by 



(19) 


where N\s the number of measurements in the sample. In essence, (19) asserts 
that increasing the number of measurements reduces the randomness of 
independently measured values, and the error will eventually converge to zero 
with an infinite number of measurements. We can express any value x being 
measured with a random uncertainty as 


X = 


x±Sx 


Ran . 


( 20 ) 


Conversely, systematic errors cannot be resolved through increased 
sampling. This is because every measurement is biased in one direction from the 
true value, and affects the value x and not the degree of uncertainty associated 
with it. Systematic uncertainties are accounted for by the addition of a correction 
factor. Correction factors may range from simple coefficients to complex higher- 
order variables, denoted by Sx^^^. However, the integration of systematic 

uncertainty into (20) is somewhat complicated, and can have three forms: 

^ = ^ + K±^Ran’ ( 21 ) 

x = x±(Sxs^^+Sx^^^), ( 22 ) 

or 

x = x±p,„f +{S^„f. ( 23 ) 

Determining which of the above equations to use is no trivial matter, and It 
should also be noted there is no rigorous mathematical proof or justification in 
selection. However, with careful consideration to each source of error; (21), (22), 
and (23) can provide a reasonable estimation into the total uncertainty of our 
measurements. Systematic errors are inherently difficult to Identify, since they 
require comparison to another (preferably more accurate) instrument. If one can 
conclusively determine than an anomalous systematic error is present with a 


28 



known value, a correction factor addition such as in (21) is an appropriate 
method to reduce the error. However, to be positively Identified, a systematic 
error’s magnitude must be significant enough to overcome other sources of error. 
As a result, there is another systematic uncertainty in the estimation of the 
correction factor. For this reason the systematic error in (21) is better 
represented by, 5x^y^=CF ±5X(,p, where the correction factor (CF) of known 

value is differentiated from the uncertainties with determining that value. 


Equations (22) and (23) can be rewritten as 


x = x + CF±(8xsy, + 5xp^^), 

(24) 

and 


A: = J + CF + ^(4,f. 

(25) 


The decision to use (24) or (25) largely depends on whether the error 
sources are independent or correlated. A simple sum of uncertainties, as in (24), 
represents the largest possible error associated with the measurement and 
would occur if these errors are 100% correlated. In many cases however, this 
may be a gross overestimation of uncertainty. If, for example, two sources of 
error are said to be no greater than 1%, the total uncertainty should be 
cumulative by use of (24) and equal 2%. However, if the given errors are said to 
have come from a normal distribution (i.e., 68% of all instruments have 
systematic errors better than or equal to 1%), and the errors are independent 
from each other, than (25) is more appropriate. In practice, exactly normal error 
distributions and complete Independence between different error sources usually 
do not exist; however, we can assume these requirements are met closely 
enough to use the standard statistical techniques for computing a reasonably 
accurate estimates of the total measurement uncertainty. The quadratic sum in 
(25) reduces the margin of uncertainty provided by (24), but also acts to properly 
weight known sources of error. In some cases, sources of error will prove to be 
insignificant when compared with other sources in (24). 
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B. INSTRUMENT ERRORS 


Each instrument has an associated accuracy and precision that 
accompany its measurements. The manufacturer usuaiiy defines these 
properties as a resuit of rigorous testing in the iaboratory. Accuracy is a 
description of how ciose a measurement is to the true value or reference value 
(systematic uncertainty). Precision, or sometimes known as reproducibility, is the 
degree to which muitipie measurements differ from each other (random 
uncertainty). Totai uncertainty is the combination of both accuracy and precision, 
and if defined by the manufacturer, may also include other sources of error in a 
variety of environmentai conditions. In the case of both the radiosonde and tower 
sensors, the manufacturers provide totai uncertainty and are iisted in Chapter III. 

According to Vaisala’s radiosonde datasheet, the totai accuracy inciudes 
repeatability, long-term stability, response time, effects due to measurement 
conditions (e.g., soiar contamination), and effects from eiectronics. However, 
Vaisaia aiso indicates that the temperature sensor performance depends on 
aititude. The response time is considerabiy faster 0.4 seconds at the surface 
compared to 2.5 seconds at 10 hPa. Repeatability is consistent at 0.15 C 
regardless of altitude, but reproducibility is much more dependable at the surface 
0.2 C, where above 10 hPa reproducibility increases to 0.5 C. Vaisaia 
engineered the radiosonde to operate within the entire verticai structure of the 
atmosphere, so ciaimed toierances must account for all conditions. Our 
experiment was limited to five meters above the ground, so the expected 
instrument error is probabiy much less than the total error provided by Vaisaia. 
Using the repeatability and reproducibility values for surface measurements an 
alternate instrument error of the radiosonde can be calcuiated 

56^^, = Vo.2"+0.15" = 0.25. (26) 

Another limitation of an instrument is the significant figures, or detaii, used 
to read measurements. Both radiosonde and tower recorded measurements from 
aii instruments to the hundredths place. Therefore, any measurement provided in 
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the experiment may be off by as much as ±0.005 units (e.g., a typical relative 
humidity recorded by tower sensor four might read 42.33%, and thus could 
actually be as high or low as 42.334% or 41.335%). The measurement of the 
seven tower sensor heights is also subject to error. These were originally 
reported in inches to the nearest y4 inch by the assembly team. We can therefore 
expect tower heights recorded to have an accuracy of 0.125 inches. When 
converted to metric this translates into ± 0.318 cm. 

C. MEASUREMENT ERRORS 

Uncertainties introduced as a result of instrument interpretation and 
experimental methods are considered measurement errors. There are two 
primary sources of measurement uncertainties in our experiment, altitude 
uncertainties while the aircraft is hovering, and contributions from the effects of 
an operating quad-copter (e.g., prop-wash and parcel mixing). 

1. Altitude Errors 

Crucial to the experiment was hovering the InstantEye at a particular 
altitude, for enough time to adequately sample the atmosphere. However, 
achieving this proved challenging for several reasons. The aircraft maintained its 
position partially through Its automatic control system, but also through operator 
control. Wind gusts created instability and altitude oscillations, as the control 
system tried to compensate. Altitude placement was also approximated by visual 
reference from the pilot since an accurate method of measurement would have 
interfered with the experiment. As a result, it is estimated that during the 
experiment, uncertainty in sUAS altitude could be as much as ±0.25 meters. 

Our altitude error is significant because it directly affects the atmospheric 

values sampled. For instance, the temperature profile is a function of height and 

varies logarithmically In an unstable atmosphere. Therefore, an altitude Increase 

of 0.25 meters will result in a potential temperature change not commensurate 

with an equal decrease in altitude from the same point. In Figure 11 is an 

example of how altitude uncertainty can affect the expected temperature 
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measurements. The red circles represent the altitude and mean potential 
temperature profile of the seven tower sensors during a flight that took place on 
May 6th. The dotted black line connecting the sensors is an interpolation of the 
measurements in log-space. The black “X” marks a distance 0.25 above and 
below each tower sensor located on the temperature profile, and the colored 
boxes are a visual representation of how much potential temperature uncertainty 
exists for the given change in altitude. It is easily observed that potential 
temperature uncertainty increases notably as altitude decreases in this very 
unstable surface layer situation. 


Vertical Temperature Profile 6 May 



Figure 11. Mean temperature profile with altitude error boxes 


Note that in Figure 11, tower level one and two are labeled as a single 
tower sensor. This is again related to altitude uncertainty. The sensors at tower 
level one and two are 11 and 19 cm above the ground, respectively. This is well 
within the 25 cm estimated altitude error, and as a result it is nearly impossible to 
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differentiate a hover between the two levels. For this reason, levels one and two 
were consolidated for the subsequent analysis. Because the temperature profile 
is approximately linear between the lower two tower sensors, we can assume 
that the average height of sensor one and two (14 cm) will also correspond to the 
average temperature recorded by both instruments. For the remainder of this 
text, measurements indicated by tower sensor one refers to the average of both 
the lowest two tower sensors. Hence, tower sensor two’s label is omitted from all 
figures and tables. 

Up to this point, we have discussed two different types of uncertainty 
related to altitude in the experiment: uncertainty in measuring the tower heights, 
and uncertainty in the aircraft altitude. Total error with respect to altitude (z) can 
be expressed as 

^z = ^(0.25f+(3.18x10"")' = 0.25m , (27) 

where the first term is the radiosonde/sUAS height error, and the second term is 
the tower elevation height error. It is clear that the error contributed by the 
uncertainty in measuring tower height is negligible when compared to the 
uncertainty of the aircraft altitude. 

Next, we must translate our error in altitude uncertainty into a useful 
parameter such as an uncertainty in potential temperature {56). We will assume 
that our altitude uncertainty is random and normally distributed, since we were 
not purposely trying to deviate from the prescribed altitude level. In Figure 11, 
height is a function of temperature, and inside the altitude error boxes we can 
say that potential temperature is approximately linear. Averaging the two 
temperatures corresponding to 5z, we get an expected value for Temperature 
with corresponding uncertainty 

6 ±56. (28) 

However, this average temperature does not correspond to the temperature 
located at the tower sensor. Instead, the expected value measured by the 
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radiosonde is slightly above or below the average tower temperature, depending 
on the curvature or second derivative of the temperature profile. Likewise, 
magnitude and uncertainty are dependent on height. This means that each flight 
has unique uncertainties that need to be calculated for each tower level. 

The magnitude of temperature error that is considered systematic is best 
represented as a correction factor while the random uncertainties are 

represented as We are assuming that our altitude estimation is random, but 

we cannot say that all measurements taken are completely Independent from 
one another; therefore, the N variable in (19) needs to be reduced from the value 
of the number of samples. This is because the aircraft takes time to recover from 
an altitude perturbation. If a wind gust brings the aircraft upward 0.25 m from a 
hovering state, it takes a few seconds for the operator and control circuitry to 
respond. During this time the sensor may make several measurements while 
elevated. The series of measurements made are related to a single altitude 
event; hence, they are not independent. We can say, however, that the group of 
measurements associated with the altitude event is independent from group 
associated with the previous and next altitude event. The application of (19) 
without grouping our samples underestimates the true random uncertainty 
associated with altitude. We know that the radiosonde records a sample every 
two seconds, and it is estimated that typical altitude perturbations require 
approximately three to five seconds to recover from. Therefore, each altitude 
event has a group of two measurements associated with it. Equation (19) can be 
rewritten for this particular example as 



Although this method is quite subjective, it is a reasonable estimate of altitude 
groups, and provides a better value of uncertainty than If we had not grouped 
them at all. 
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2 . 


Aircraft Errors 


A second source of measurement uncertainty is the InstantEye platform. 
The four rotors spinning above the meteorological sensors likely have some 
effect on the surrounding air, but the extent to which this influences sensor 
measurement was not known before this research was performed. A major 
objective of this thesis work is to provide an answer to this question. This can 
only be achieved by accounting for all other sources of error, and incorporating 
the total uncertainty to the final analysis. 

D. ENVIRONMENTAL ERRORS 

The experiment was conducted outside in a non-controlled environment. 
As a result environmental influences further contributed to the total uncertainty in 
our measurements. Specifically, solar contamination and natural turbulent 
fluctuations were significant enough to be identified in the data. 

1. Solar Contamination 

Solar contamination is a significant problem with meteorological sensors. 
Direct exposure to the sun causes warming of the sensor’s components, and in 
turn can affect the quality of measurements taken. The tower sensors were 
carefully designed to minimize solar influence through a double wall sensor 
enclosure and the installation of a solar shield, so It is assumed to be negligible. 
Conversely, both temperature and relative humidity are of particular concern on 
the RS92 since their sensing components are positioned on a metallic boom 
protruding outward from the main chassis directly exposed to the sun. Relative 
humidity is affected more by solar contamination due to its dependence on 
temperature. Several noteworthy attempts have been made to accurately 
account for solar contamination with the application of a relative humidity 
correction for the RS92 based on solar angle (Cady-Pereira 2008, Vomel 2007, 
Yoneyama 2008). While this correction is integrated into the RS92 software 
making it transparent to the user, the precision of the relative humidity sensor is 
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still inadequate for the conditions observed in our experiment, the details of 
which are further addressed in Chapter V. 

The temperature probe is to a lesser degree affected by solar 
contamination, but still significant enough to contribute to the total uncertainty of 
temperature measurements. According to Vaisala (2010), the temperature 
sensor correction at sea level can be as much as 0.1 C at the maximum solar 
angle. This is also compensated for in the RS92 software, but requires two key 
assumptions in the operation of the radiosonde; the sensor boom is positioned at 
an angle, and it is attached to a weather balloon. Both of these measures ensure 
there is adequate ventilation over the sensor components. 

A typical weather balloon ascends vertically at a velocity of about 5 m/s, 
and the air passing over the boom minimizes solar contamination effects. In our 
experiment, the entire radiosonde was positioned horizontally from its normal 
(balloon) operating position, negating the requirement for the boom to be angled. 
When the InstantEye was hovering there was no vertical motion of the sensor; 
however, eddies produced by the rotors assisted in the ventilation requirement. 
In the experiment, we were able maintain the InstantEye at ground level for 
approximately a minute with the rotors operating. We then turned the rotors off 
for an additional minute and observed an average temperature rise of 3.1 C. This 
is clear evidence that the propellers serve a key function in reducing solar 
contamination by providing ventilation to the RS92 sensors. 

We then examined temperature observations with the radiosonde placed 
in the direct sunlight and compared them to those made in the shade. In both 
cases the InstantEye was hovering to ensure the rotors were providing sensor 
ventilation. In both observations no discernable difference in temperature 
measurements were observed, suggesting that ventilation is adequate over the 
sensor boom to reduce solar contamination to the manufacture’s engineering 
specifications. This is confirmed in other studies conducted by Jauhiainen et al. 
(2014) and Moti (2014), where temperature comparison between day and night 

balloon soundings revealed a reproducibility of less than 0.1 C, with a standard 
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deviation of 0.02 C at the surface. We can safely conclude that the uncertainties 
contributed by solar contamination, while attached to an operating InstantEye, 
are sufficiently negligible that they do not need to be added to the manufacturer’s 
claimed accuracy, which incorporates this error into its measurement algorithms. 

2. Turbulent Fluctuations 

In an unstable atmosphere, it is expected to observe energy transfer from 
rising air parcels in the surface layer as a result of the ground being warmer than 
the atmosphere. These motions create localized turbulence and mixing, which 
can introduce random errors into our measurements of mean temperature 
values. Fortunately, these short-term variations can be reduced through 
sampling; however, one must first be able to isolate the contributions due to 
turbulence from all the other sources of noise. For this reason, it is quite difficult 
to develop a strategy to reduce uncertainty in our final measurements. 

To get a feel for the data collected during the experiment, a time series 
plot is presented (Figure 12) over a fifteen-minute period for a flight conducted in 
May. This particular flight took place in the morning after sunrise in the presence 
of an unstable temperature profile. The reference tower measurements are 
displayed in multicolored lines with each color representing a unique sensor on 
the tower. Tower level one in dark blue is closest to the ground, and tower seven 
in red is the highest sensor on the tower. The tower sensors captured a series of 
peaks and troughs in potential temperature taking place over several minutes. 
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7 May 2014, Starting at 1637Z 



Figure 12. Time series plot of temperature with large variations 

Potential temperature captured by the radiosonde while in flight (black 
line) is overlaid on top of the data from the various tower sensors. This particular 
15-mlnute section of flight was chosen as an example since the aircraft was in 
operation over the entire period, and was attempting to hover at altitudes 
corresponding to tower measurement levels. The variation in potential 
temperature observed by the radiosonde Is quite extreme in both magnitude and 
time scale, and uncorrelated on small time scales with the variations captured by 
the tower. This can partially be explained by the fact that aircraft altitude and 
associated temperature error is absent from Figure 12. Another explanation can 
be related to the proximity of the tower and aircraft to each other. During a flight 
the aircraft was flown at a distance of around five to 10 meters away from the 
tower. This served to prevent a collision between the aircraft and tower, and also 
to prevent contamination of the tower sensors due to any prop-wash effects 
generated by the aircraft. Any turbulent eddies smaller than the separation 
distance produced uncorrelated temperature fluctuations. 
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7 May 2014, Starting at 1637Z 



70 75 80 85 

Time (min) 

Figure 13. Time series plot of potential temperature with aircraft heights 

Figure 13 is the same time series as Figure 12 but with altitude data 
incorporated. The solid black line still represents the radiosonde’s potential 
temperature record, but with the addition of color-coded data points. Each point 
is displayed with a unique shape and represents a hovering altitude. The 
color associated with the data is same color as the corresponding tower level 
(e.g., when the aircraft was said to be hovering at level five, all the recorded 
measurements are indicated by a yellow square—yellow being the color 
representing tower five). 

There are portions of the radiosonde’s profile In Figure 13 without 
accompanying altitude data. These reflect periods of altitude transition, or times 
when control of the aircraft was difficult. These data were rejected in the analysis 
since it is unclear where the aircraft was located at these times. However, when 
looking at the data of value (I.e., those points containing altitude information) the 
variability is still very large for a given hovering level. Since the tower sensors 
had less variability than the aircraft, there are three possible explanations: 
(1) This is a result of altitude variation, (2) the variations are natural In origin but 
are not adequately being captured by the tower sensors due to a slower 
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response time, or (3) the forces and motions of the aircraft rotors are driving 
large-scale variations. 

The temperature sensors on both the tower and radiosonde have some 
key differences in sampling methods. The radiosonde uses a fast response 
capacitance sensor, which records a temperature reading every two seconds. 
The tower sensors use an RTD, which average over a 30-second time interval, 
recording temperature values every 10 seconds. By design, the tower 
temperature sensor will smooth out large short-term temperature fluctuations if 
present. To determine if the same variations recorded by the radiosonde were 
also present at the tower, a fast response thermocouple was installed at tower 
level four during the August and November JIFX flights. In doing so, the 
thermocouple is devoid of any altitude uncertainty or adverse effects contributed 
by the aircraft rotors. The thermocouple also records its temperature 
measurement every second, providing an improved sample rate over both the 
tower sensors and the radiosonde. 

Figures 14 and 15 provide sample windows comparing the three different 
types of temperature sensors under two different conditions. The unstable 
atmosphere is presented in Figure 14, taking place during a very warm August 
afternoon. The stable atmosphere is shown in Figure 15, which took place on a 
cool November morning before sunrise. Note, due to a processing error, 
radiosonde data were only recorded to within 0.1 C for the stable profile. In both 
cases, large short-term variation in potential temperature is captured by both the 
thermocouple and radiosonde. The magnitude of variation observed in the 
unstable atmosphere is larger when compared to the stable atmosphere. More 
importantly, the radiosonde and thermocouple recorded temperature variations 
nearly identical in magnitude for both cases. Since there was no altitude or prop¬ 
wash uncertainty in the thermocouple measurements, the same fluctuations 
captured by the radiosonde must also be natural in origin. This is because any 
additional variations on the radiosonde measurements induced by other effects 
such as prop wash would be expected to cause higher variability. Other sources 
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of error included in the radiosonde’s measurements are not easily identifiable in 
Figures 14 and 15. Figure 16 is a spectral plot of both the thermocouple and 
radiosonde showing similar results, including the -2/3 slope that is associated 
with turbulence in the inertial subrange. The apparent increase in spectral 
magnitude at the highest frequencies is due to aliasing; no filtering was 
performed. This similarity is further evidence that the radiosonde is observing 
fluctuations natural in origin. We would expect more high frequency variations if 
they were originating from the prop-wash. 


a) Tower level 4 RTD (6 samples/min) 



22 23 24 25 26 


d) c) Aircraft hovering at levels 3, 4, and 5 (30 samples/min) 



Time (min) 

Figure 14. Time series of potential temperature (unstable conditions) 
taken on August 12th, 2014 starting at 2126Z, 
a) tower level four RTD, b) tower level four TC, and 
c) aircraft hovering at levels three, four, and five. 
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Potential Temperature (6) 
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c) Aircraft hovering at levels 3, 4, and 5 (30 samples/m in) 



Time (min) 


a) Tower Level 4 RTD (6 samples/min) 



Figure 15. Time series of potential temperature (stable conditions) taken 

on November 5th, 2014 starting at 1413Z, 
a) tower level four RTD, b) tower level four TC, and 
c) aircraft hovering at levels three, four, and five. 
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Figure 16. Spectral plot of potential temperature August 12th, 
thermocouple (top) verses radiosonde (bottom). 


The presence of natural temperature variations aiso reaffirms our choice 
in using the seven slow-response RTD tower sensors for anaiysis. At any instant 
in time, the measured temperature at the thermocouple will not likely match the 
temperature measured by the radiosonde as a resuit of turbulent fluctuations. 
The difference could be off by as much as four degrees, as in the exampie of 
Figure 14. In this case, a time-running average, as provided inherently by the 
tower RTD sensors is best to describe the mean state of the atmosphere in the 
presence of this natural variability. The seven sensors on the tower in effect 
perform this function through their longer response time averaging out short-term 
variations. This provides an accurate reference mean temperature for 
comparison with the radiosonde. 
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Since turbulent fluctuations are pseudo-random, their uncertainty can be 
reduced through sampling. However, determining the variation of fluctuations is 
particularly challenging in this experiment due to the absence of data. The 
measurements recorded by the radiosonde while in flight capture turbulent 
fluctuations as shown above, but also contain altitude and prop-wash 
uncertainty embedded within. To determine the variability of natural fluctuations 
an independent fast-response probe is required, and our thermocouple was 
installed late in the experiment. There are only three flights with fast-response 
data available, and four flights without. 

Immediately following one flight in August we used the radiosonde to 
gather some turbulent variation data. During this test, the InstantEye was shut 
down and manually held next to the lowest five tower sensors for approximately 
one minute. This removed any unwanted altitude, prop-wash, or spatial 
uncertainty in the measurements when compared to each tower sensor. Tower 
levels six and seven were not included in the test because they were too high to 
reach and keep stable for consistent measurements. Table 4 summarizes those 
results. 


Level 

N 

<l 

b 

5 

30 

0.58 

4 

30 

1.03 

3 

30 

1.15 

1 

47 

1.21 


Table 4. Difference in potential temperature test results 
(radiosonde-tower) with props off, 
and held manually. 


The results of the test demonstrate that standard deviation decreases with 
height and is in alignment with surface layer theory. According to an empirically 
derived study (Garratt 1992), temperature variations in the surface layer during 

unstable conditions can be represented as a function of height (z), 
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where 6'*is the temperature scaling parameter, and I is the Monin-Obukhov 
length scale. If we assume that turbulent fluctuations are consistent, then 6'*and 
L are by definition constant in the surface layer, so (30) can be represented as 

(^ffrurb=cz\ (31) 

where c is a constant for a given 0* and L. Equation (31) shows how variation 
is expected to diminish with height, and our test results confirm a close 
relationship. Therefore, not only do we require an independent fast-response 
sensor during each flight, one is also required at each level to accurately 
determine the variation of natural fluctuations. We fall well short of this 
requirement so must rely on analytic methods to account for these uncertainties. 

We should also expand upon our independent test conducted at the 
lowest five tower levels to include the fast-response thermocouple. Looking at the 
difference between the thermocouple and tower level four for the entire flight, our 
sample size is increased to 6650 and the standard deviation of potential 
temperature difference between tower and thermocouple is 1.25 C. Since our 
test with the radiosonde was limited to 30 samples producing a standard 
deviation at level four of 1.03 C, it is clear that we need a larger sample size to 
adequately represent the variation of turbulent fluctuations. Knowing the heights 
of each tower level, and having sampled 0„, =1.25 C from the thermocouple, 

^^Turb 

we can then solve for c in (27) and apply the equation to other levels. This helps 
resolve our turbulent fluctuation uncertainty at different heights, but does not 
account for the flights where the thermocouple was unavailable for a reference. 

A larger unstable temperature gradient will produce greater temperature 
variability (Garratt 1992). Given this relationship, we can estimate the standard 
deviation of short-term temperature fluctuations by only knowing the difference in 
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temperature. This is possible since the tower was left to record measurements 
continuously for several days, and accumulated an abundance of data from the 
thermocouple and tower sensor four. Figure 17 shows the relationship between 
temperature gradient and variation. Over four days, the thermocouple was 
compared to the potential temperature recorded by tower sensor four. Standard 
deviation was calculated every 15 minutes using an hour window (3600 
samples), and then plotted against the average temperature gradient recorded 
during that hour (tower one-tower seven). The results clearly show a linear 
relationship between the two variables, which can be represented by 

<7„„^=0.40(g,-g,) + 0.05. (32) 



Figure 17. Relationship between temperature gradient and standard 
deviation of the difference in potential temperature between the 
thermocouple and tower sensor 4 


Using (28) we can generate a reasonable estimate of cr., , the standard 

^^Turb 

deviation due to turbulent fluctuations at level four. Equation (27) can then be 
used to estimate the variation at each level knowing only the average 
temperature gradient for unstable atmospheres. For stable atmospheres, we do 
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not expect to see variations as a function of height. Garratt confirms that a 
relationship between a-g and other parameters under stable conditions are 

uncertain. Under stable conditions we know that small variations appear in the 
data as shown in Figure 15, so we can assume similar variations were also 
present during the May 8th morning flight. For these reasons, we would be safe 
to assume similar variations at each tower level. Examining the uncertainty from 
stable conditions in the August and November data, 95% of cr., samples from 

the thermocouple remained below 0.37 C. This figure is probably larger than 
what actually occurred that morning, but gives us a safe value of uncertainty for 
analysis. 

Similarly to our altitude uncertainty, we cannot assume that 
measurements with turbulent fluctuations are independent from each other. Each 
fluctuation requires time to pass through a parcel. We estimated that this time 
was about six seconds, or three measurements from the radiosonde. For 
turbulent fluctuations. Equation (19) can be rewritten as 


^0Turb 



(33) 


Tables 5 and 6 present the potential temperature gradients, uncertainties, and 
uncertainties in the mean used for calculation of errors associated with turbulent 
fluctuations. 


Level 

May 6 

May 7 

May 7 

May 8 

Aug 12 

Aug 13 

Nov 5 

7 

17.86 

14.00 

19.50 

5.90 

31.28 

22.33 

5.05 

1 

21.28 

17.46 

23.20 

5.59 

34.46 

25.21 

3.41 

A0 

3.42 

3.46 

3.70 

-0.31 

3.18 

2.88 

-1.64 


Table 5. Average flight temperatures in C recorded 
at levels one and seven 
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May 6 

May 7 

May 7 

12 Aug 

12 Aug 


A0 

=3.42 

Ad 

=3.46 

Ad 

=3.70 

Ad 

=3.18 

Ad 

=2.88 

Level 

^Turb 

^Turb 

^Turb 

^Turb 

^Turb 

^Turb 

^Turb 

^Turb 

^Turb 

^Turb 

7 

0.73 

0.13 

0.74 

0.17 

0.79 

0.11 

0.68 

0.10 

0.62 

0.15 ■ 

6 

0.93 

0.13 

0.94 

0.15 

1.00 

0.16 

0.86 

0.14 

0.79 

0.16 ; 

5 

1.14 

0.17 

1.16 

0.18 

1.23 

0.24 

1.07 

0.15 

0.97 

0.16 1 

4 

1.43 

0.33 

1.45 

0.15 

1.55 

0.22 

1.34 

0.16 

1.22 

0.16 

3 

1.79 

0.47 

1.81 

0.24 

1.93 

0.38 1 

1.67 

0.22 

1.52 

0.36 i 

1 

I 2.40 

1.86 

2.43 

0.35 

2.59 

0.39 

2.24 

0.52 

2.04 

0.35 : 



May 8 

Nov 5 


A6» = 

-0.31 

Ad = 

-1.64 

Level 

^Turb 

^Turb 

^^Turb 

^Turb 

7 

0.37 

0.12 

0.26 

0.05 

6 

0.37 

0.11 

0.26 

0.04 

5 

0.37 

0.09 

0.26 

0.03 : 

4 

0.37 

0.08 

0.26 

0.04 

3 

0.37 

0.12 

0.26 

0.04 

1 

0.37 

0.12 

0.26 

0.04 


Table 6. Estimated uncertainty in potentiai temperature variation 
and estimated uncertainty in the mean j 
due to turbulent fluctuations. 


E. BIAS ERRORS 

Some of our data revealed evidence of the radiosonde being significantly 
biased as it was recorded during flight. In these cases, we chose to discard our 
results because of the large uncertainty. However, we only conducted two flights 
under stable conditions and both of these contained a large degree of bias. 
Because it is important to anaiyze both stable and unstable atmospheres we 
needed to correct for this uncertainty in the best way possible. 


48 






















Figure 18. Stable time series November 5th uncorrected 



Figure 19. Stable time series November 5th corrected for bias errors 

Figures 18 and 19 provide an example of a bias error recorded on the 
morning of November 5th just before sunrise. The radiosonde temperature 
measurements are about 2.4 C cooler than the tower record as shown in Figure 
18. The slope over the time series indicates that a gradual warming event 
is underway, and its being captured by both instruments suggests a simple 
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correction factor is an appropriate solution to the discrepancy. To accomplish 
this, we assumed temperatures recorded at level seven were most accurate due 
to a weaker gradient. The data were then fit to reflect temperatures recorded at 
level seven similar to those recorded by the tower instrument, as seen in red on 
Figure 19. 

There are several possibilities why this bias was introduced. It could be 
that the radiosonde was not calibrated sufficiently. During the experiment 
radiosondes were used multiple times, and as a result could be subject to 
calibration issues. Another possibility could be the result of physical disturbance 
to the radiosonde. Our flights required some skill at keeping the aircraft airborne 
especially at lower levels, and the temperature sensor brushed against the 
ground during several “crashes” and hard landings over the course of the 
experiments. The data confirmed that on occasion a new bias error would occur 
after one of these events. It is uncertain if vibrations changed the internal 
calibration, or if the sensor was chemically or materially contaminated by 
accumulated debris in the form of soil or moisture. 

F. ERROR SUMMARY 

After examining the various categories of uncertainty in our 
measurements, we can now formulate an acceptable method to interpret our 
measurements. We will focus on potential temperature since it is the primary 
parameter we are evaluating. Since we are comparing the temperature recorded 
by the radiosonde to the tower’s measurement, we will normalize our 
measurements and perform statistical studies on the difference {^0). Through 
averaging we can determine by 

^e = T9+CF,,^^ + CF,,±^e (34) 

where is the correction factor due to bias errors, is a correction factor 

due temperature offset related to altitude uncertainty. Tolerance is given by 50 
defined as 
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(35) 


SO56\^^ + + 56\^^ + 56].^^^ ^ 

where and SOj.^^ are the instrument errors associated with the radiosonde 
and tower, respectively, 50^^^ is the uncertainty in temperature due to altitude, 
and is the uncertainty related to turbulent fluctuations, all assumed to be 
independent. Equations (34) and (35) are then applied to the data set, 
formulating Table 7, providing total correction factors and uncertainty for each 
flight at each tower level. 


Level 

May 6 

^^Total SO 

May 7 

^^Total SO 

May 7 

^^Total SO 

12 Aug 

^^Total SO 

12 Aug 

CFiotal SO 

7 

-0.01 

0.41 

0.00 

0.43 

-0.01 

0.41 

0.00 

0.40 

0.00 

0.42 , 

6 

0.00 

0.41 

0.00 

0.42 

0.00 

0.42 

0.16 

0.42 

-0.01 

0.42 1 

5 

-0.02 

0.43 

-0.04 

0.43 

-0.01 

0.46 

0.03 

0.42 

0.01 

0.42 

4 

-0.08 

0.52 

-0.10 

0.42 

-0.14 

0.45 

-0.34 

0.43 

-0.05 

0.42 : 

3 

-0.39 

0.64 

-0.30 

0.47 

-0.56 

0.58 

-0.33 

0.47 

-0.37 

0.56 

1 

0.39 

1.96 

0.51 

0.53 

0.48 

0.56 

1.06 

0.71 

0.53 

0.54 


Level 

May 8 

^^Total SO 

Nov 5 

^^Total 36 

7 

-0.75 

0.41 

2.59 

0.39 : 

6 

-0.92 

0.41 

2.34 

0.39 

5 

-0.72 

0.40 

2.35 

0.39 1 

4 

-0.90 

0.40 

2.41 

0.39 : 

3 

-0.90 

0.41 

2.40 

0.40 

1 

-0.77 

0.41 

2.21 

0.39 


Table 7. Total correction factor and uncertainty 
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V. ANALYSIS AND RESULTS 


A. HUMIDITY 

Humidity is an important atmospheric property affecting the characteristics 
of EM propagation, and our anaiysis wouid not be compiete without addressing 
this in our findings in the experiment. Unfortunateiy, due to the dry surface at 
Camp Roberts, specific humidity gradients (cioseiy reiated to the vapor pressure, 
which controis EM refraction) were negiigibie during each of our fiights, so 
performing the same type of anaiysis that was performed for potentiai 
temperature was not possibie. Tabie 8 provides the mean specific humidity in 
g/kg recorded by the bottom six sensors for each of the fiights. Reiative humidity 
was the actuai property measured by the sensors, but since it is dependent on 
temperature, reiative humidity was converted to specific humidity, a conserved 
quantity described in Chapter ii. The dataiogger did not record humidity 
measurements from tower ievei seven so these data are absent. By examining 
the quantities depicted in Tabie 7, no obvious trends exist as a function of height. 
Differences between tower ieveis were aiso insignificant, and the maximum 
difference observed between any two tower ieveis was no greater than 0.50 g/kg 
for each of the fiights. 


Level 

May 

6 

May 

7 

May 

7 

May 

8 

Aug 

12 

Aug 

13 

Aug 

13 

Nov 

4 

Nov 

4 

Nov 

4 

Nov 

5 

Nov 

5 

6 

4.45 

5.66 

5.68 

4.99 

4.45 

5.66 

5.68 

6.35 

5.90 

5.09 

5.13 

4.97 

5 

4.66 

5.77 

5.88 

4.99 

4.66 

5.77 

5.88 

6.40 

6.00 

5.28 

5.13 

4.94 

4 

4.63 

5.74 

5.76 

4.89 

4.63 

5.74 

5.76 

6.46 

6.04 

5.31 

5.23 

4.98 

3 

4.54 

5.80 

5.65 

4.97 

4.54 

5.80 

5.65 

6.58 

6.12 

5.20 

5.20 

4.90 

1 

4.67 

5.98 

5.85 

4.95 

4.67 

5.98 

5.85 

6.76 

6.40 

5.33 

5.18 

4.83 

Maximum 

Difference 

0.22 

0.32 

0.23 

0.10 

0.22 

0.32 

0.23 

0.41 

0.50 

0.24 

0.09 

0.16 


Table 8. Mean values for specific humidity (g/kg) by tower ievei 

and flight 
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A result of relatively homogeneous moisture content is the inability of the 
radiosonde to differentiate between measurements made at different altitudes. 
Figure 20 depicts the relationship between relative humidity as a function of 
temperature and specific humidity in a standard atmosphere of 1000 hPa using 
typical values observed during the experiment. The accuracy of the radiosonde’s 
relative humidity sensor is ±5%. Figure 20 demonstrates that a 5% uncertainty 
in relative humidity corresponds to an uncertainty in specific humidity, around 
0.5 g/kg. Since our observed differences in absolute humidity were mostly less 
than 0.5 g/kg we cannot analyze humidity given the accuracy of the radiosonde 
and lack of substantial gradients. For these reasons, we will focus our analysis 
on temperature only. Because potential temperature and specific humidity are 
both scalar quantities that are conserved and not affected by the sUAS except by 
air movement, we believe that the sUAS effects on sampling (air movement from 
different levels and turbulence) Is the same for specific humidity as for potential 
temperature. 
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B. POTENTIAL TEMPERATURE 


We conducted a total of 12 flights during the three JIFX. Of these, we 
rejected five flights from our analysis due to obvious temperature bias errors or 
lack of sufficient data points. From the remaining seven flights, five were 
conducted in unstabie conditions while the sun was present, and two were 
conducted in stabie conditions before sunrise. The two flights conducted in a 
stable atmosphere also contained bias errors, but due to the absence of other 
reliable data, correction factors were used to account for errors as described in 
Chapter IV. 

Our temperature analysis was oniy possibie because we were able to 
compare measurements observed from the InstantEye to those observed from 
the tower at a corresponding aititude. We can mathematicaiiy represent our 
comparison by ietting s represent the ievei at which the radiosonde was 
hovering, and t represent the tower ievei for each measurement /. The mean 
difference of potentiai temperature {AO) then can be expressed as 

= where 5,rG{l,3,4,5,6,7}, (36) 

^ ’ i=\ 

where N\s the total number of measurements observed during any single flight. 
When s and fare equai the expected vaiue of AO^^ shouid equal zero if no errors, 

environmentai variations or sensor differences are present. 

1. Unstable Profile 

Resuits conducted from the five fiights during unstable temperature 
conditions are shown in Table 9. In this particular table, 5 and t equal means 
that the radiosonde measurements are oniy being compared to the tower ievei in 
which the InstantEye was hovering adjacent to. The number of measurements is 
represented by Af, is the mean difference in potential temperature, is 

the standard deviation of difference in potential temperature, and 56 is the 
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expected uncertainty of our calculation of A6' derived in Chapter IV. Note, that 
there were two separate flights that took place on May 7th. 




a) A0 

6 May 


Level 

N 

AO 

*^A6' 

SO 

7 

94 

0.43 

0.83 

0.60 

6 

155 

0.50 

0.97 

0.60 

5 

141 

-0.10 

0.84 

0.61 

4 

55 

0.47 

1.08 

0.68 

3 

44 

0.16 

1.05 

0.77 

1 

5 

0.08 

0.44 

2.01 



b) AO 

7 May 


Level 

N 

AO 

*^A6' 

SO 

7 

56 

0.09 

0.69 

0.61 

6 

119 

0.23 

0.76 

0.60 

5 

121 

0.25 

1.00 

0.61 

4 

288 

0.36 

1.26 

0.60 

3 

175 

0.15 

1.44 

0.64 

1 

147 

0.17 

1.15 

0.69 



<1 

7 May 


Level 

N 

AO 


SO 

7 

158 

-0.26 

0.65 

0.59 

6 

125 

-0.23 

0.86 

0.60 

5 

81 

-0.01 

1.15 

0.63 

4 

150 

-0.19 

1.11 

0.62 

3 

76 

0.22 

1.36 

0.72 

1 

135 

0.44 

1.36 

0.71 


d) 12 August 


Level 

N 

AO 

Q5 

<1 

b 

SO 

7 

139 

0.03 

0.82 

0.59 

6 

115 

-0.03 

0.82 

0.60 

5 

158 

0.10 

0.87 

0.60 

4 

202 

-0.27 

1.28 

0.61 

3 

167 

-0.44 

1.59 

0.64 

1 

56 

0.25 

1.29 

0.83 


e) A^ 13 August 

Level 

N 

AO 

Q5 

<1 

b 

SO 

7 

52 

0.17 

0.49 

0.60 

6 

69 

-0.28 

0.80 

0.61 

5 

112 

-0.11 

1.02 

0.61 

4 

180 

0.18 

1.08 

0.60 

3 

53 

-0.07 

1.01 

0.71 

1 

103 

■ 34l 

0.98 

1 069 1 


Table 9. Difference in potential temperature statistics of the 
radiosonde at hover compared to equivalent tower level in 
unstable conditions, a) May 6th, b) May 7th, c) May 7th, 
d) August 12th, and e) August 13th 
56 


There are some general trends observed from Table 9. First, each 
calculated mean A0 is well within its corresponding derived tolerances values 
S0 with the exception of August 13th at level one (shaded). This suggests that 
the InstantEye’s effect on its radiosonde payload measurements was not 
detectable among the other sources of uncertainty and error. However, this 
tolerance is suspected to be high, since S0 was calculated using the total 
instrument error claimed by the factory, ±0.5 C. Recalculation of S0 with the 
alternate instrument error (26) derived in Chapter IV Is shown in Table 10. 


Level 

6 May 

A0 56 

7 May 

56 

7 May 

A6 56 

12 Aug 

A6 56 

13 Aug 

A6 56 

7 

1 0 43 

O. 4 H 

0.09 

0.43 

-0.26 

0.41 

0.03 

0.40 

0.17 

0.42 

6 

1 0 50 

0.4 J 

0.23 

0.42 

-0.23 

0.42 

-0.03 

0.42 

-0.28 

0.42 

5 

-0.10 

0.43 

0.25 

0.43 

-0.01 

0.46 

0.10 

0.42 

-0.11 

0.42 

4 

0.47 

0.52 

0.36 

0.42 

-0.19 

0.45 

-0.27 

0.43 

0.18 

0.42 

3 

0.16 

0.64 

0.15 

0.47 

0.22 

0.58 

-0.44 

0.47 

-0.07 

0.56 

1 

0.08 

1.96 

0.17 

0.53 

0.44 

0.56 

0.25 

0.71 

|T34^ 



Table 10. Mean difference in potential temperature between the 
radiosonde at hover compared to equivalent tower level in 
unstable conditions with alternate radiosonde instrument error 


With a smaller tolerance for uncertainty, levels six and seven are now 
suspect measurements for the May 6th flight, but this is not a reoccurring trend 
for other flights. Most of the mean potential temperature values still remain within 
our expected tolerance. Our results still suggest that InstantEye contamination of 
potential temperature measurements Is undetectable given our equipment and 
environment. 

Another general trend observed from Table 9 is a decrease in the 
radiosonde temperature measurement variation (cr^^) with respect to height. This 

means that the accuracy of our measurements is greatly improved at the highest 

tower levels. This can be explained by the decreasing temperature gradient at 
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higher levels, as well as the related reduction in turbulent fluctuations as the 
aircraft moves away from the surface. The combination of these two factors also 
reduces our expected tolerance from the error Equation (35) as altitude 
uncertainty and fluctuations become increasingly insignificant above five meters, 
leaving instrument errors as the primary source of measurement uncertainty. If 
only instrument errors are considered, (35) would equal 0.39 C, and Table 10 
indicates that 50 \s converging to a number very close this value at tower level 
seven (4.93 m above the surface). 

We suspect that the prop-wash effects would cause the sampled air 
parcels to come from a different level than the actually location of the sUAS. To 
quantify this scenario. Figure 21 is presented showing the mean radiosonde 
measurements compared to every other tower level. For instance, the mean 
potential temperature observed by the radiosonde while hovering at level four, 
can be compared with all six mean tower observations. The diagonal from bottom 
left to upper right are the figures 50 from Table 10, 5 and t being equal. Figure 
21 is color-coded corresponding to the magnitude of mean potential temperature 
difference, and the column containing the value closest to zero is circled. Since 
the radiosonde measurements were subtracted from the tower’s measurements, 
a positive value indicates the radiosonde is on average measuring warmer 
temperatures than the tower it is being compared to. Likewise, negative values 
indicate the radiosonde is measuring on average cooler temperatures than the 
tower sensor it is being compared to. 

When examining each column of Figure 21, the 50 value closest to zero 
is an important location since it corresponds to the tower level and altitude. The 
zero point within the column implies that the mean air being sampled is 
originating from this altitude. One would expect to observe the mean difference 
closest to zero when 5 and t are equal (i.e., along the diagonal) if there is no 
influence from the InstantEye on vertical displacement of the sampled air. The 
circles in Figure 21 correspond to the tower levels closest to the zero value. This 

implies that circles above the diagonal are sampling air being pulled down from 
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the rotors, originating from a higher altitude. Conversely, circles below the 
diagonal imply that the mean air sampled is being pulled up from below as a 
result of rotor motions. The latter could occur if air is recirculated upward around 
the downwash. Circles along the diagonal suggest that the average air sampled 
is not originating from elsewhere, but is consistent with the atmosphere at its 
hovering altitude. 


a) A6/ 6 May 



1 3 4 5 6 7 

Aircraft Level (s) 


b) AO 7 May 



1 3 4 5 6 7 

Aircraft Level (s) 


d) A4^12 August 



1 3 4 5 6 7 


Aircraft Level (s) 


e) A0^3 August 



1 3 4 5 6 7 


Aircraft Level (s) 


c) AO 7 May 



1 3 4 5 6 7 

Aircraft Level (s) 


Figure 21. Mean potential temperature difference, radiosonde (s) 
compared to each tower sensor (t). 
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Figure 21 illustrates that the potential temperatures from the tower and 
radiosonde match quite well along the diagonal, meaning the air being sampled 
has not been mechanically mixed sufficiently to significantly affect the mean 
values. Table 11 summarizes the positions of the circles corresponding to closest 
zero value in Figure 21. Each level is tallied into three categories; above, below, 
or on target (i.e., on the diagonal), totaling five for each of the separate flights. 
With the exception of level six, most levels indicate a strong tendency for the 
radiosonde to sample air similar to its corresponding tower level. Where this is 
not true, differences in /s.0 are slight, and a pattern cannot be established. There 
is no evidence at any altitude that the radiosonde tends to sample air above or 
below the rotor blades when averaged over a period of time. This is somewhat 
surprising as one might expect that the sampled air would always come from 
above the sUAS. It appears that some air that was pushed down by the props is 
recirculated back up outside of the downwash zone and brings some lower air up 
with it. 


Level 

1 

3 

4 

5 

6 

7 

Above 

1 

1 

1 

1 

2 

— 

On Target 

4 

4 

3 

4 

0 

3 

Below 

- 

0 

1 

0 

3 

2 


Table 11. Tally of mean potential temperature difference above, below, 

or on target with respective tower level. 


2. Stable Profile 

There were only two flights conducted in the presence of a stable 
temperature profile. For this reason, the data for these flights were not rejected 
due to bias errors prevalent in the time series, but instead correction factors were 
applied. This creates difficulties in determining where the sample air comes from, 
since consistent sampling of air above or below would also create biases in the 
temperature. But we can look at the overall characteristics of the comparisons, 
even if the correction factor may not be exact. Results for the two flights are 
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summarized in Table 12, with the measurement uncertainty {56) calculated with 
the smaller radiosonde Instrument error. 


a) 59 8 May b) 59 5 November 


Level 

N 

59 

<1 

b 

59 

Level 


59 

<1 

b 

59 

7 

31 

0.05 

0.12 

0.41 

7 

97 

0.00 

0.09 

0.39 

6 

34 

0.25 

0.25 

0.41 

6 

106 


0.15 


5 

52 

0.01 

0.13 

0.40 

5 

216 

0.03 

0.25 

0.39 

4 

64 

0.05 

0.06 

0.40 

4 

132 

0.10 

0.12 

0.39 

3 

27 

0.18 

0.12 

0.41 

3 

108 

0.20 

0.19 

0.40 

1 

28 

0.15 

0.06 

0.41 

1 

131 

Lzl 

0.39 

1 0 39[ 


Table 12. Difference in potential temperature statistics of the 
radiosonde at hover compared to equivalent tower level in 
stable conditions, a) May 8th, b) November 5th 


The differences in potential temperature are within the expected 
uncertainties, the exceptions being level one and six on November 5th. Note that 
the sample size from the flight on May 8th is much smaller than November’s, due 
to limited flight time before sunrise. A smaller sample size may be insufficient to 
determine accurate means. 

Unlike the unstable profile, there is no evidence that uncertainty in our 
measurements (cr^^) decreases with height. Also the expected variation of our 

measurements {59) remains approximately constant in a stable profile. This is 
explained by reduced turbulent fluctuations that do not vary with height 
(discussed in Chapter IV) and a smaller temperature gradient resulting in less 
altitude-temperature uncertainty. The error tolerance is primarily from instrument 
errors, which will not reduce further. Similar to the unstable profile analysis 
above, a comparison of the radiosonde’s mean potential temperature to each 
tower level is presented in Figure 22. 
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Figure 22. Mean potential temperature difference between the 
radiosonde hovering at a particular level compared to each tower 
sensor for two separate flights in stable conditions, a) May 8th, 

b) November 5th 

Again, the circles in Figure 22 represent the minimum absolute difference 
in potential temperature, and their arrangement is not along the diagonal as they 
were In Figure 21. Furthermore, choosing a value of ls£ closest to zero is almost 
futile, because there exist several values very close (±0.1 C) to zero within a 
column. Most of the values in Figure 22 also fall within their expected 
uncertainties, approximately 0.40 C. 

As expected, the difference in potential temperature at tower level seven 
compared to the radiosonde hovering at seven is near zero. This is because the 
bias correction factors were formulated to ensure this was the case. There is 
clear evidence from Figure 19 in Chapter IV that a distinct temperature gradient 
is observed by the radiosonde between levels seven and five, but not at the 
lower levels. This can again be verified by column seven of Figure 22, indicating 
a distinction between observations taken at level seven compared to air sampled 
by the tower below it. 

Our results also indicate that, below level seven, the radiosonde appeared 
to sample a layer of uniform temperature, even though we know a gradient 
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actually was present. Figures 19 and 22 illustrate the difficulty in differentiating 
the measurements taken while the radiosonde was hovering at levels five, four, 
three, or one. This can be explained by the InstantEye’s prop-wash forcing the 
atmosphere to mix below level seven. Since the atmosphere is stable, and the 
lower levels are cooler than the upper layers, a major source of heat energy is 
absent from the surface to restore equilibrium to the local parcel of air the 
InstantEye is operating in. In this scenario, advection from adjacent air carried by 
the wind is the dominant force to restore equilibrium. However, it appears that the 
wind was too light (<3kts) during the morning flights to negate the mixing effects 
of the InstantEye. 

When a stable temperature layer is mixed the cooler air on the surface will 
become warmer, and the warmer air aloft will become cooler; Table 12 shows 
this happening. In both cases equals zero around level five. Below level five, 
gradually increases, meaning the mean temperature observed by the 
radiosonde is warmer than mean temperature observed by the lower tower 
levels. Figure 23 depicts the mean potential temperature profile for the flight on 
November 5th. While the aircraft is hovering, the props force aloft air downward, 
displacing and mixing with the ambient and lower air. If the air being forced down 
confronts a boundary, in this case the surface, displacement and mixing are 
strengthened both laterally and upward since downward motion is limited and 
mass must be conserved. In the bottom five levels we observe a relatively 
homogeneous atmosphere, air is being forced from aloft near level six and 
mixing the lower levels. This explains the unusually high observed mean 
potential temperature at level one in Table 12. 

The temperature gradient between levels five and six are approximately 

neutral in Figure 23. This prevents a drastic change in A0 while hovering at level 

five, since the air aloft is of similar temperature. However, the temperature 

gradient between levels six and seven is steeper. At some altitude above level 

six there is an inflection point and steepening temperature gradient, though it is 

not apparent from our data where it is. The surface is no longer as restrictive to 
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flow as it was while hovering at lower levels, so mixing is predominantly occurring 
between levels five, six, and the warmer layer aloft. This explains the 0.46 C 
warmer than expected A0 at level six in Table 12, suggesting the air is 
originating as high as 2 meters above the InstantEye. 



Potential Temperature (^) 


Figure 23. Mean potential temperature height profile (first 15 minutes) 
observed from the tower sensors November 5th in 
(stable conditions) 

3. Turbulent Fluctuations 

Our contrasting results between temperature measurements taken in 
stable and unstable atmospheres require a revisit to temperature fluctuations. 
We observed that in the presence of strong turbulent fluctuations, the prop-wash 
was on average an insignificant factor for temperature measurement accuracy. 
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However, in the presence of weak turbulent fluctuations, average temperature 
measurements indicate the action of the rotors are introducing contamination by 
mixing the vertical air column. Table 13 is a comparison of the standard deviation 
of the mean potential temperature observed with the expected standard deviation 
contributed solely by turbulent fluctuations. Greyed boxes indicate the observed 
variations in temperature measurements exceed those expected from turbulent 
fluctuations. 



Level 

May 8 

^ffTurb 

Nov 5 

^ffTurb 

7 

0.12 

0.37 

0.09 

0.26 i 

6 

0.25 

0.37 

0.15 

0.26 

5 

0.13 

0.37 

0.25 

0.26 i 

4 

0.06 

0.37 

0.12 

0.26 

3 

0.12 

0.37 

0.19 

0.26 i 

1 

0.06 

0.37 

0 39 



Table 13. Standard deviation of mean potential temperature observed 
{(T^g), and expected standard deviation due to turbulent 

fluctuations 
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In the presence of an unstable atmosphere we expect to see variations in 
temperature due to turbulent fluctuations; strong near the surface and decreasing 
with height. This relationship is observed in our measurements. However, most 
of the variation in our observations is less than what was expected, suggesting 
that the InstantEye is reducing variations caused by turbulent fluctuations. It is 
only at levels five, six, and seven that we observe the variation of our 
measurements exceeding those expected from turbulent fluctuations, thus other 
forces are contributing to uncertainty at higher altitudes. Also note that as our 
observations approach the surface variations in measurements become larger, 
but not to extent we expect from the growing contributions due to turbulent 
fluctuations. Furthermore, the growing uncertainty of our measurements 
produces an anomalous reversal at level one where is usually less than level 

three. 

For the flights in a stable temperature profile, uncertainty from turbulent 
fluctuations has no height dependency, so we assume they remain constant with 
elevation. The uncertainties in our measurements are smaller than what we 
observed for the unstable case due to smaller magnitude gradients and turbulent 
fluctuations, but an unusually large uncertainty is observed at level one for the 
November flight. A similar finding did not appear in the May flight, but this may be 
related to the limited number of samples taken (28 in May versus 121 in 
November). This result shows the InstantEye making a strong contribution to 
temperature measurement variations near the surface. 

Our findings in potential temperature observations indicate that the 
InstantEye contributes to measurement uncertainty while hovering, and this 
contribution is greatest near the ground. Turbulent fluctuations were present in 
every flight, but the InstantEye’s contamination assisted in counteracting this 
source of uncertainty. 

As the InstantEye gains altitude, variations due to turbulent fluctuations 
decrease in an unstable atmosphere, and eventually the InstantEye’s 


66 



contributions to measurement uncertainty become apparent. This explains why 
on August 13, the radiosonde measured more variation than was expected due 
to fluctuations at level five (Table 13). The temperature gradient on this day was 
the least of the five unstable flights, and in turn variations from turbulent 
fluctuations were also the least. On this day the InstantEye’s effect on 
measurement uncertainty was observed at a lower level than usual. Note how on 
both May 7th flights, the temperature gradient was the largest, and turbulent 
fluctuation variations were also the largest of the seven unstable flights. On these 
flights, fluctuations were so strong that the InstantEye could not overpower them 
at level seven. 

As the InstantEye descends to the layer just above the surface (level one) 
prop-wash contamination effects become stronger, as do turbulent fluctuations. 
In the stable case fluctuations were comparatively small, so variations in the 
measurements hovering at level one was dominated by InstantEye 
contamination. Whereas in the unstable flights turbulent fluctuations were 
strongest at level one, but the amplified effect from the InstantEye helped to 
reduce measurement variations considerably. In the case of August 13th the 
observed variation In temperature measurements was smaller than the other 
flights at level one. The reduced affect from natural fluctuations together with the 
amplified affect from the InstantEye established a close to equilibrium state. 
Evidence of this can be found in the anomalous rise of potential temperature 
difference (Figure 21) at level one. The surface’s energy transfer is restricted 
causing the air in the vicinity of the InstantEye to warm. 

C. PRESSURE 

Another important atmospheric parameter for meteorological and EM 
prediction is pressure. The pressure differential across a rotor plane generates lift 
(Layton 1984), and therefore we should expect to observe higher pressure from 
the radiosonde since it is positioned beneath the rotor plane. To quantify the 
magnitude of pressure error due to lift, we examined measurements taken during 
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the August 13th flight. Near the end of the flight we hovered for a period of time 
at the lowest six levels. Shortly after, we turned the aircraft off and held the 
radiosonde at the lowest six tower levels for a period of time. Level seven was 
too high to reach by hand. Figure 24 and Table 14 illustrate the data collected 
and the results. 
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Figure 24. Dynamic pressure captured by the radiosonde on August 
13th, a) InstantEye in flight, b) InstantEye not in flight 
(held at each level by hand). 



I In Flight 


Held by hand 


Level 

P 

ap 

P 

ap 

AP 

6 

982.23 

0.11 

982.05 

0.04 

0.18 

5 

982.38 

0.08 

982.17 

0.07 

0.21 

4 

982.38 

0.08 

982.25 

0.07 

0.13 

3 

982.42 

0.08 

982.30 

0.05 

0.12 

1 

982.59 

0.12 

982.36 

0.06 

0.23 


Table 14. Mean, standard deviation, and mean difference of pressure 
recorded by the radiosonde August 13th, with and without 

rotors operating. 
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The data points in Figure 24 exhibit a degree of uncertainty for both cases. 
When the radiosonde was held by hand, there is no altitude uncertainty, and 
since we are only comparing results against the same instrument, we only need 
the radiosonde’s instrument error to determine the expected uncertainty. Vaisala 
claims that at the surface repeatability of the pressure sensor is within 0.4 hPa 
assuming a 95% confidence interval (i.e., ±2cr where (t= 0.1 hPa), and the data 
presented in Figure 24b as well as Table 14 indicate this agrees with the present 
observations. 

The radiosonde in flight must account for altitude uncertainty in addition to 
instrument error. We estimated our altitude error to be ±0.25 m, which translates 
into an expected pressure uncertainty of 0.03 hPa, since the expected rate of 
pressure change with elevation in a standard dry atmosphere is approximately 
0.12 hPa/m (Petty 2008). Total pressure uncertainty from the radiosonde in flight 
becomes 

SP = ylo.f+0.03^ =0.104 hPa. (37) 

This is a negligible difference, and while in flight, all but the pressure 
measurements from levels six and one remained within expected uncertainty. 

The differences between mean pressures values (AP) in Table 14 
suggest that there is a constant pressure bias from the InstantEye ranging from 
0.12-0.23 hPa. This corresponds to a 1-2 m change in altitude, and there is no 
observable relationship between pressure bias and tower level. In a similar 
experiment. Dr. Peter Guest of the Naval Postgraduate School found pressure 
bias to increase near the surface (Figure 25) while in flight due to ground-effect. 
This effect became apparent at about 1.3 m, and the maximum mean bias was 
observed to be 0.6 hPa just above the surface near 16 cm (Guest and Machado 
2015). 

Altitude oscillations during our flight may explain our inability to recreate 
Figure 25. Above we showed that expected pressure change due to altitude 
variations are negligible with an uncertainty of ±0.25 m, but we did not take into 
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account the InstantEye’s effect on pressure beneath the rotors. The aircraft uses 
its internal pressure sensor to maintain altitude. If the aircraft starts to descend, 
the internal pressure sensor will detect an increase in pressure and signals the 
InstantEye to increase thrust. This creates an artificial pressure increase beneath 
the rotors lifting the aircraft back to the desired altitude. This was a common 
occurrence during our experiment, and usually resulted in overshooting the 
desired altitude, oscillating up and down until equilibrium was reached. This 
phenomenon was exacerbated by wind gusts, turbulence, and operating near the 
ground in the presence of ground effect, but a time-running average helped 
mitigate this in our temperature analysis. Since our sample size is small, we 
would expect to see altitude oscillations appear in the results. Large oscillations 
would translate into more variation in the pressure measurements and produce a 
less accurate mean. Observing data points of Figure 25 suggests that while the 
InstantEye was in flight very few altitude oscillations were occurring since the 
variation is similar for both cases. This is not the case if Figure 24, suggesting 
altitude oscillations were significant. 



Minutes Minutes 

Figure 25. InstantEye pressure effect from Guest’s experiment (shown 
negative as a proxy for elevation), left-1 nstantEye in flight, 
right-held by hand props off (courtesy of Dr. Peter Guest) 


Our results do present an interesting finding: the mean pressure increase 
by the radiosonde in flight was at least 0.12 hPa above the mean pressure when 
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held by hand at each level. Data in Figure 25 at level six appear to confirm this 
pressure bias when in flight. This means while hovering at any altitude we can 
expect to observe a 0.12 hPa pressure measurement bias, translating into a 1 m 
altitude bias. Guest’s observations suggest that pressure will further increase 
below 1.3 m as the surface is approached due to ground-effect. While we did not 
observe this to be true in our observations, large altitude oscillations may explain 
this inconsistency. 
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VI. CONCLUSIONS 


A. SUMMARY OF EXPERIMENT 

This experiment is one of the first studies to analyze quad-rotor effects on 
meteorological sensors. Our use of an RS92 radiosonde onboard the InstantEye 
sUAS produced some important results. In unstable temperature profiles, 
temperature measurements made In the surface layer are as accurate as the 
manufacturers claimed accuracy for surface measurements. In stable 
temperature profiles mixing occurs below 1.3 m, smoothing out the temperature 
gradient making it difficult to differentiate between levels. Above 2 m in stable 
conditions, temperature measurements were warmer than expected; suggesting 
air is originating from as high as 2 m aloft. The InstantEye’s rotors contribute to 
the variation in temperature measurements, and this effect is strongest near the 
surface. The variations introduced by the prop-wash help suppress variations by 
turbulent fluctuations, but are present above the surface layer. The InstantEye 
also introduces a 0.12 hPa pressure bias while In flight due to the pressure 
effects of the rotors corresponding to a 1 m elevation bias. While we could not 
quantify the presence of ground-effect in our results, other experiments have 
shown ground-effect to introduce a pressure bias of 0.60 hPa (5 m altitude error) 
near the surface. 

B. OPERATIONAL APPLICABILITY 

This experiment must also be considered in the correct context. Our flights 
were conducted in extreme meteorological conditions, since we were confined to 
the surface layer in a desert-like environment. The atmosphere was dry, and 
strong temperature gradients dominated lower levels during unstable flights. As a 
result, turbulent fluctuations were intense during unstable temperature 
conditions. In a real naval application, our InstantEye-radiosonde combination 
would be used at sea, where there is obviously a surface humidity source, and 
temperature gradients are generally not as strong. As a result turbulent 
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fluctuations of temperature would be less prominent, but humidity fluctuations 
would be stronger. 

Our experiment also required hovering the InstantEye for periods of time 
in one location. This was to obtain sufficient data for statistical analysis since our 
reference tower was limited to the lowest five meters. Taking multiple 
measurements while hovering helps reduce uncertainty and Improve accuracy, 
but It is suspected that hovering also gives time for rotor contamination to take 
hold mixing the ambient air. In practice, the aircraft and sensor would ascend to 
an altitude of about 600 to 900 meters to sample a large vertical profile in the 
boundary layer. Constantly moving upward is advantageous for meteorological 
sensing because It also reduces the Induced velocity through the rotors, as we 
saw in Chapter II. 

C. RECOMMENDATIONS FOR FUTURE RESEARCH 

There were several limitations to our experiment that warrant further 
investigation. Trials above the surface layer are a logical next step forward to 
investigate the effects of prop-wash on vertical motion, both upward and 
downward, and measurement uncertainty in the absence of strong turbulent 
fluctuations. Our results suggest that measurements taken in an unstable 
atmosphere will perform very well above the surface layer, but it is unclear if rotor 
mixing will significantly bias temperature and humidity accuracy in stable 
atmospheres, or if there is a dependency on the strength of the gradient. 

Sea trials are necessary to observe the effects of prop-wash on humidity 
since we did not have a strong humidity gradient in our experiment. Flights 
originating from ships would also reveal operating considerations and limitations 
for future use in a naval weather-sensing program. Furthermore, other 
meteorological sensors could be used In place of Vaisala’s RS92, as well as an 
alternative to the InstantEye. There are numerous combinations of quad-rotors 
and sensors that could be used in place of our equipment to optimize sensor 
placement and accuracy, as well as identify an ideal cost-benefit relationship. 
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